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ABSTRACT 



We consider a correlated wavefunction including particle-hole pairing at half 
a reciprocal lattice vector for itinerant electrons hopping on a square lattice in 
two dimensions and subject both to on-site and nearest-neighbor repulsion. We 
find, within mean field theory, that in a suitable range of parameters there is a ten- 
dency toward a novel form of ordering characterized by discrete symmetry breaking 
and a gap at the magnetic zone boundary, which supports hole-like semiconductor 
behavior near half filling despite an apparently normal Fermi surface. The effec- 
tive quasiparticle couplings to phonons includes a significant d-wave piece, which 
plausibly leads to an especially robust tendency toward superconductivity. 
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One has become accustomed, especially in connection with the superfluid 
phases of He 3 , to the existence of quite intricate pairing correlations between par- 
ticles. On the other hand, well-known forms of ordering including ferromagnetic, 
antiferromagnetic, and charge-density wave can be described in terms of particle- 
hole pairing [1,2]. It is quite natural, then, to consider by analogy the possibility 
of states of matter characterized by particle-hole pairing of more general types. In 
this note we shall show that indeed one such ordering is suggested to be favor- 
able in an analysis of the mean-field correlation energy induced by very simple, 
purely repulsive interactions. This may have significant implications for the Mott 
insulator problem, as we shall discuss. 

1. Analysis 

Let us consider electrons on a square lattice in two dimensions with spacing a. 
We suppose that the screened Coulomb interaction can be modelled by an on-site 
repulsion of magnitude U and nearest-neighbor repulsion of magnitude V, that is: 

Hint. = U^c\(x) C] {x)c\(x) Ci {x) + ^cS(i)^(i)cJ(a/)^(a/), (1.1) 

x x,x' 

where x' runs over the four nearest neighbors of x, and a, (3 =j, j are summed 
over. Transforming (1.1) to momentum space, we find the form 

Hint. = U ^2c\(ki)c^(k 2 )c\(ks)ci(h) 8 L {k\ + h - k 2 - k 4 ) 

+ ci{ki)c a {k 2 ) c^(k 3 )cp(k 4 ) C(k 3 - k 4 ) 5 L (h + k 3 -k 2 - h) , 

hi 

(1.2) 

where 5i denotes the delta function in the reciprocal lattice, and C(k3 — k 4 ) = 
cos(/c3 — k^) x a + cos(/c3 — k^) y a. 

Interesting possibilities for particle-hole correlations, whose energetic conse- 
quences can be seen in a mean- field analysis of (1.2), involve ordering at momen- 
tum transfer Q, where 2Q is a reciprocal lattice vector. We shall focus on the case 
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Q = {ji/a,-n/a). In the spirit of BCS theory, we postulate a wave-function of the 
form 

where the product runs over some subset of the interior of the magnetic zone (i.e. 
the diamond (k x ± k y )a = ±7r), and \u\ 2 + \v | 2 = 1 for all values of the indices, to 
preserve normalization, and the number of factors is determined by the density of 
electrons (see below). The expectation value of H mt . in the state \I/ is 

(tf|i?int.|*) = ^2 U ( U kl V kl +U k ]V k] ){uk'[V k '[ +Uk'[Vk'[) 
k,k' 

- 2V{u k] v k] + MfcT^fcT + Uk[V k [ + u k[ v ki ) 

X (Mfc'T^fe'T + u k']V k ^ + U k 'lV k 'l + Uk'lVk'l) 

+ VC(k-k') ((Mfct^T - «fcT^t)(^'T^'T -^fc'T^fc't) 

+ (uklVkl ~ UklVkl){uk'lVk'l - Uk'lVk'l)) , 

(1.4) 

where the bar denotes complex conjugate. 

Without loss of generality we may take the us to be real and non-negative, and 
write v ka = w ka exp(i(f) ka ) with the ws real and non-negative. Then (1.4) becomes 

(*|i?int.|*) = ^2 4f/ ( M fcT w fcT M fe'I w fe'|) COS0 fc |COS0 fc /| 
k,k' 

- 8V(u k ^W k ^ COS0fc| + U k [W k [ COS0fc|)(Mfc/|Wfc/| COS0fc'| + U k 'lW k 'l COS0fc/|) 

- AV C(k - k') ((M fcT w fcT sin fc T M fc'T w fc'T sin <t>k"\) 

+ ( u klWklsin(f) kl u k/l w k/i sm(f) k/i )) . 

(1.5) 

At this point, it is easy to identify three essentially different ways whereby 
favorable correlation energy can emerge: 

• From the first term, which favors | cos0| = 1 with opposite signs for spin up 
and spin down - this corresponds to antiferromagnetism (this is the case analyzed 
in [1]). 
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• From the second term, which favors | cos0| = 1 with the same sign for spin 
up and spin down - this corresponds to a charge density wave. 

• From the third term, which favors |sin0| = 1 (that is, v pure imaginary 
relative to u). 

The third possibility is in many ways the most interesting, and we shall now 
discuss it in some detail before returning to compare it with the other two. 

To continue the analysis we must minimize the full energy, including both one- 
particle (hopping) and correlation parts. For the one-particle part, we take the 
nearest-neighbor form e k = — t(cos k x a + cos k y a). For simplicity let us consider 
nearly half filling, so that the chemical potential is small and the Fermi surface 
is near the boundary of the magnetic zone. The equations for the different spin 
components do not affect one another, so we shall temporarily drop the spin index. 
By standard methods one then finds the solution of the minimization problem in 
the form 




(1.6) 



where e k = \{e k - e k+ q) and 



A fc = AV ^2 u k'W k ' sin fc sin 4> k > C(k - k') 
k> 



(1.7) 



satisfies the gap equation 



sin <p k sin <p k / C(k — k') . 



(1.8) 




and C(k - k' 



) can be written in the more revealing form 
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The energy will be minimized when A is non-negative; together with the form 
of the gap equation this determines A k to be of the form A^ = A\ sin k x a\ (or, what 
is equivalent energetically, A k = B\ sin k y a\) which has a p- wave-like symmetry on 
the Fermi surface or A k = C\ cos k x a\ (= C\ cos k y a\ on the diamond) which is 
d- wave-like. The sign of sin0^ will adjust itself, changing at the zeroes of A so 
as to enforce these forms. The gap equation is quite complicated in general, but 
simplifies near half filling. One has in this case 



TV 





for the A form and 



l = 4V I^Lf di - Sin I „ (1.10) 



7T 



i = %v[*L[ di - cos I . (i.n) 

J sine J J~e2 + C 2 cosH 

Q V 

The factor, which reflects the van Hove singularity in the density of states 
exactly at half filling, favors the C form, as does the factor 2 arising from the 
equality of the first two terms in the final factor of (1.8). Away from half-filling, 
or in the presence of a next-nearest neighbor hopping term, there will not be a van 
Hove singularity at the Fermi surface. However, an enhanced density of states will 
remain if these effects are not too large. At half filling one has occupation up to 
the edge of the diamond, where e vanishes, and thus a solution to the gap equation 
- even apart from the van Hove enhancement - for arbitrarily weak coupling. 

To discuss the situation below half filling, we must consider how the condition 
on electron density is implemented. The operators 

At =- t _- t [ ' 

°ka — V kaC ka U ka C k+Qa , 

analogous to the familiar Bogoliubov-Valatin operators in superconductivity the- 
ory, are a complete set of fermion creation and annihilation operators with diagonal 
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anticommutation relations. 7; and 5; create quasiparticles below and above the 
gap, respectively, or, in semiconductor terminology, in the the valence and conduc- 
tion bands. At half-filling, all of the 7^ , or conduction band, states are filled and 
the ground state (1.3) can be written \1> = Yl k T k ^T k ^\0 >, where k ranges over all 
momenta in the magnetic zone. In general k ranges only over momenta within the 
Fermi surface (actually a curve, in two dimensions), which is determined by the 
condition + A| = /i. The chemical potential, /i, is chosen so that the area en- 
closed by the Fermi surface is equal to the desired density. Below half-filling the de 
integral will be cut off by /i instead of diverging as C — > 0. Thus the gap equation 
will not have a non-trivial solution at arbitrarily weak coupling, but only starting 
at a finite value, dependent on the deviation from half-filling, of the coupling. 

Nothing in our considerations so far has correlated spin up and spin down, 
so that spin singlet {u k ] = u k \), spin triplet (u k ^ = — lifcj,), or intermediate pos- 
sibilities are equally favorable. The degeneracy will be split if we include into 
our Hamiltonian a nearest-neighbor spin-spin interaction Js ■ s where of course 
s(x) = c a (x)aoC^(x) is the electron spin operator. When such a term is included 
in the Hamiltonian, contributions sensitive to the pairing correlations arise from 
the crossed channel. A short calculation shows that for J > 0, antiferromagnetic 
coupling, the spin singlet state is lowered in energy, and the triplet raised. Of 
course the assumed effective coupling here might have either sign, and it need not 
be at all the same as the effective J used in a different approximate description of 
the same material (e.g. in at — J model). In any case the spin singlet order will be 
favored by its relative immunity from fluctuations, as we shall presently discuss. 

Now let us compare the energies of the competing orders. These are quite com- 
plicated and depend on the interplay between the single particle energies and the 
interaction terms (1.4) (with the u k s and v k s appropriate to these orders substi- 
tuted). Some qualitative features of the phase diagram may be obtained from the 
coefficients of the terms in the purely repulsive Hamiltonian H int . We find that the 
terms favoring the three types of states (antiferromagnetic, charge density, "C-type 
d-density") are in the ratio —U : U — 8V : — 4V(cos 2 k x a). Here the last factor 
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arises from the non-trivial angular dependence of the d-density ordering and the 
pairing interactions that give rise to it. It will be substantially larger than ^ for the 
C form of d-density ordering, because the density of states is largest where the co- 
sine is unity (indeed, to the extent the van Hove singularity dominates, this factor 
is unity), though its precise value is highly model- and doping-dependent. When 
this factor is unity, the C form of d-density ordering is stable only at U — AV, 
according to this naive criterion. However, one should also consider that antifer- 
romagnetic ordering, since it breaks a continuous symmetry, is subject to severe 
fluctuation effects - in principle, for example, a two-dimensional system is rigor- 
ously forbidden to exhibit such ordering at any non-zero temperature - which can 
severely degrade the favorable mean-field energy. In a realistic domain of parame- 
ters for CuC>2 planes Schrieffer et al. [1] found that the magnetic moment density 
is renormalized by a factor .6 due to fluctuations, which is also approximately the 
factor indicated by experiment. One might expect the correlation energy to be 
renormalized roughly by the square of this factor. The singlet d-density order, as 
we shall discuss further below, breaks only a discrete symmetry, so it is safer in this 
regard. Also the d-density order, because it contains nodes, has a more favorable 
one-particle energy. The one-particle energy for this state, in the presence of a 
gap of magnitude A, is 6.40A 2 /t while the antiferromagnet or charge- density wave 
have one-particle energy 10.65A 2 /t. (These numerical coefficients were obtained 
by integrating the one particle energies up to a chemical potential /i = 2A which 
avoids the van Hove singularity. The results are similar for other choices of the 
chemical potential; the robust point is that the d-density is substantially favored.) 
Altogether, then, it seems that d-density order, especially in its singlet form, is a 
serious candidate to describe real states of matter. 



8 



2. Comments 



1. Since neither the Hamiltonian we analyzed nor the approximation scheme 
we employed is exact, it behooves us to identify what feature of the proposed state 
might be expected to have precise validity. As in many cases in physics, a broken 
symmetry is the heart of the matter. One finds in the spin-singlet d-density state 
the non-vanishing expectation value 

(4(fc + Q)cP{k)) = i5?f(k) , (2.1) 

where / is a real function that changes sign under a n/2 rotation. The % indicates 
that time-reversal symmetry T is violated, and the angular dependence indicates 
that symmetry under 7r/2 rotation is violated. Furthermore, the dependence on 
the momentum offset Q = (7r/a,ir/a) indicates that symmetry under translation 
through a single lattice spacing is also violated, since this too changes the sign 
of the order parameter. Let us swiftly add, however, that one can combine any 
two of these broken symmetry operations to recover a valid symmetry. Thus, in 
particular, one does not expect to see any simple direct macroscopic manifestation 
of the symmetry breaking. 

Indeed, it is not entirely easy to identify accessible signatures for d-density 
order even microscopically. There are no extra coherent peaks in elastic x-ray 
or neutron scattering, since the relevant single-particle expectation values cancel 
upon integrating over angles. Perhaps the most fundamental signature derives from 
the existence of low energy particle-hole excitations at momenta connecting points 
on the Fermi surface repeated after translation through Q. These should reveal 
themselves in inelastic neutron scattering, as an anomalous extended continuum. 
Closely related structure would also be expected in photoemission. 

Although the model we have analyzed above is two-dimensional, the essential 
idea of the ordering pattern (2.1) is not intrinsically two-dimensional, and invites 
generalization to three-dimensional materials. 
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2. We have already sketched the construction of the quasiparticles. When 
one replaces the normal electron operators by these quasiparticles, what emerges 
near half filling can be described to a first approximation using a semiconductor- 
like picture. A gap has opened up at the boundary of the magnetic zone, so 
that doping below half-filling it will produce an effective concentration of hole-like 
charge carriers proportional to the doping. The Fermi surface will have its normal 
area or volume, though if the nominal Fermi surface passes through the magnetic 
zone boundary it will be distorted from its nominal band theory form. The critical 
behavior at the transition will depend crucially on whether or not this occurs [3,4]. 
Transport in the spin-triplet d-wave state will be anomalous, due to the Nambu- 
Goldstone mode associated with breaking of the spin rotation symmetry. (More 
precisely, there will be important fluctuations, but not necessarily a simple Nambu- 
Goldstone mode, in the appropriate channel.) Transport in the spin-singlet d-wave 
state will also be anomalous, due to fluctuations associated with transition to this 
state at a finite coupling and doping. The precise nature of these anomalies is 
under study. 

Because it effectively halves the size of the Brillouin zone without introduc- 
ing any manifest antiferromagnetic or charge density ordering, the d-density wave 
ordering appears to be an interesting candidate to describe Mott insulators (or 
the related doped "semiconductors") which exhibit neither antiferromagnetic nor 
charge- density ordering. 

3. Because the quasiparticles of the d-density state are composed of d-wave 
mixtures of fundamental electronic modes, the effective phonon-mediated interac- 
tion between these particles will contain a substantial d-wave component. Since 
it is less sensitive to Coulomb repulsion, effective attraction in this partial wave 
plausibly leads to a particularly robust instability toward d-wave superconductiv- 
ity. (Of course, it is still subject to the strong fluctuations associated with two- 
dimensionality for an isolated layer.) The accumulation of density of states near 
the magnetic zone, which is close to the Fermi surface for small doping, enhances 
this tendency. 
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4. In the states we have considered above, the mean density of spin up and 
spin down have been supposed equal. Actually the [/-term energy (though not 
the single-particle energy) is minimized by a completely polarized, ferromagnetic 
state. Within the polarized state one can still examine the possibility of d-density 
ordering, which we find can indeed be energetically favorable. This provides a 
model for possible ferromagnetic Mott insulator behavior at \ filling. 

Acknowledgements: We wish to thank R. Schrieffer and S. Sondhi for valuable 
comments. 
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